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Abstract 

The formalism in order to obtain the Dark Energy equation of state is 
extended to non-flat universes and we consider the inequalities that must be 
satisfied by Phantom Dark Energy in this case. We show that due to a non- 
vanishing spatial curvature satisfying the observational bounds, the uncertainty 
on the determination of the Dark Energy equation of state parameter w, when 
it is taken constant, can be significant and that it is minimal for some redshift 
z cr ~ 3. We consider the potential of future measurements of the gravitational 
waves emitted by binaries at high redshifts z > z cr to reduce this uncertainty. 
Results obtained here should also be relevant for a weakly varying equation of 
state with w ~ — 1. 
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1 Introduction 



There is growing observational evidence supporting the idea that the universe ex- 
pansion departs radically from the standard picture pQ. Assuming a spatially flat 
isotropic universe, about two thirds of the critical energy density seems to be stored 
in a component called Dark Energy, able to produce an accelerated expansion. The 
precise nature of Dark Energy is still unknown and the subject of intense investiga- 
tions and speculations The main reason is that the simplest solution still fitting 
the data, a cosmological constant A, is not satisfactory as this implies an extraordi- 
nary fine tuning of the vacuum energy. On one hand it has to be unnaturally small - 
about 123 orders of magnitude smaller than its "natural" value - and still not zero, on 
the other hand observations themselves may force us to consider Dark Energy with a 
varying equation of state jS] . A popular and extensively studied possibility is a mini- 
mally coupled, slowly rolling, scalar field called Quintessence which would reproduce 
a mechanism quite similar to that of inflation 4J. In these models, it is possible to 
obtain a negative pressure and accelerated expansion through an adequate balance of 
kinetic and potential energy. They satisfy the Null Energy Condition p+p > 0, which 
is saturated by a cosmological constant A, p\ = —p\- However, a further surprise 
may come from the observational requirement to have Phantom Dark Energy with 
p + p < 0. 

We have access to the expansion history of our universe at low redshifts z using 
Supernovae data. By measuring their luminosity distances di{z) we can reconstruct 
the Hubble diagram H(z) with a differentiation procedure. To have access to the 
equation of state parameter w(z) = ^|4, one has to go through another differentiation 
prodecure, so the data do not allow to determine w(z) sharply at the present time. 
On the other hand, this equation of state can give insight into the microscopic nature 
of Dark Energy. This is especially true if Dark Energy is of the Phantom type 
as it considerably restricts the possible candidates. Indeed, it is well known that 
the intensively studied quintessence models where the Dark Energy sector is some 
minimally coupled scalar field cannot account for an equation of state with w < — 1. 
To account for this possibility, an extension of these models was originally suggested 
with the sign of the kinetic energy opposite to the conventional one [5]. Interestingly, 
we note that scalar-tensor theories of gravity provide Dark Energy models which allow 
for Phantom Dark Energy 111,05 and where in addition the latter will be (weakly) 
clustered [8J. If observations force us to adopt Phantom Dark Energy, a whole class 
of models will be ruled out. It is further known that the future of our universe is 
dramatically different with Phantom Dark Energy i6J. Hence the special role played 
by the bordercase w = —1, also called the Phantom divide [Zj. 

It is therefore of particular interest to investigate whether the observations will 
allow us to determine formally if Dark Energy (DE) is of the Phantom type or not. 
This is bound to become a crucial issue if observations gradually select a viable region 
in the vicinity of w = — 1. So we would like to investigate the influence of priors con- 
cerning the geometry on the conclusions to be drawn from the observations. Usually, 
investigations assume a spatially flat universe. While such an assumption certainly 
has well founded theoretical motivations coming from the inflationary paradigm, it is 
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clear that a small amount of curvature is still allowed by the Cosmic Microwave Back- 
ground (CMB) data, even when combined with other data [101,1111 • As we show in 
details, the assumption of spatial flatness i.e. Qk,o exactly zero, may lead us to confuse 
Phantom DE with usual DE. This kind of cosmic degeneracy or cosmic confusion is 
something recurrent in the interpretation of data. Clearly, the full problem requires a 
comprehensive statistical analysis involving all the cosmological parameters and gen- 
eral equations of state (see e.g.[5j.[I2].[TB].[Hj). However, previous analyses assume 
spatial flatness and we want here to single out the degeneracy due to a non vanishing 
spatial curvature. As a first step, insight can be gained by considering a constant w 
[Tlj and working in the (w,Qk,o) parameter plane. We show that measurements of 
the luminosity distances at high redshifts z > 3 can substantially restrict the degen- 
eracy in the (w,flkfl) parameter plane in the neighbourhood of (w = — l,fi fc0 = 0). 
This shows the potential of high- 2: 6l measurements when combined with low- z SNIa 
data by SNAP/JDEM (Supernova Acceleration Probe/ Joint Dark Energy Mission) 
[To] and we think this is a stimulating possibility worth to explore, especially in view 
of the future space mission LISA (Laser Interferometer Space Antenna) ^0] for the 
detection of low frequency gravitational waves. 



2 The equation of state in non-flat universes 

In this Section we derive the basic equations related to the DE equation of state as 
well as the inequalities that can decide whether DE is of the Phantom type or not, 
when extended to non-flat universes. We assume our universe is isotropic with metric 

ds 2 = dt 2 - a{t) dt 2 , (1) 

and its (expansion) dynamics is fully encoded in the time evolution of the scale factor 
a(t) obeying the usual Friedmann-Lemaitre-Robertson- Walker equations of General 
Relativity. The expansion of the universe can be probed through the measurement 
of the luminosity-distance d^{ Z^j clS db function of redshift z (we take c = 1) 



d L {z) = (i + z) h, 1 \n k , \-* s (\n ki0 \* J* -^j ■ 



(2) 



Here S(u) = sinu for a closed universe, S(u) = sinhw for an open universe while S is 
the identity for a flat universe. In eq.(0), h(z) stands for the (dimensionless) reduced 
Hubble parameter h(z) = H = H(z = 0), and for small z it is given by 

h 2 ( z ) = fi m , (1 + zf + n DEfi f( z ) + n kfi (1 + zf , (3) 

, 3H 2 

where flo = f° r an Y component while f2 fc0 = — and p cr0 = 3— %. The 

pcr,0 ' &q -Hq ' c57T(_r 

unknown function f(z) expresses the evolution of Dark Energy with the expansion 

/(*) - ^ , (4) 
Pde,o 
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and is directly related to the equation of state parameter w(z) = v °^u\ ■ Indeed, the 
energy conservation equation valid for any isotropic perfect fluid 

f t =-3H(p + P ), (5) 

where t is the cosmological time, leads straightforewardly to 

, l + w{z')~ 



f(z) = exp 



dz 



1 + z' 



(6) 



For constant equation of state, we recover the well-known result f(z) = (1 + z) 3( - 1+w \ 
An important conclusion can be drawn from (JHJ), or (jSJ), for Phantom DE (w < — 1): 
its energy density decreases with increasing redshifts or increases in the course of 
time. This would give rise to a singularity in the future, the so called Big Rip. 

Though there is strong theoretical motivation from the inflationary paradigm in 
favour of a flat universe, we consider here spatially curved spaces in the conservative 
limits allowed by observations. The most stringent bound on the spatial curvature 
comes from the CMB data and these can be made even tighter in combinination with 
other data. They favour a spatially flat universe but some spatial curvature is still 
allowed. Actually the uncertainty on the spatial curvature increases with our lack of 
knowledge of DE, more specifically of its equation of state which is exactly what we 
dont know yet. 

The observational signature for Phantom DE is that the following inequality be 
satisfied 

dh 2 

— < 3 Q mfi (1 + zf + 2 fl fc>0 (1 + z) . (7) 

We see immediately from (j7J) that the assumption of spatial flatness can lead to 
an erroneous conclusion: even though (JJJ) is satisfied assuming flatness, it could no 
longer be the case if the universe is closed. Similarly, if (j2J) is not satisfied assuming 
flatness, it could still hold for an open universe. Assuming the following observational 
uncertainties 

a < a mfi <c (8) 

-B < n kfi <D, (9) 

where all constants A, B, C, D are positive, we are assured DE is of the Phantom type 
if the following inequality is satisfied 

^ <3 A [l + zf-2 B (1 + z) . (10) 

Note that (f4,o) mm = —B is negative and corresponds to a closed universe. On the 
other hand we are assured DE is not of the phantom type if 

dh 2 

— >3C (l + z) 2 + 2 D (1 + z) . (11) 
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For all other cases 



dh 2 

3 A (1 + z) 2 - 2 B (1 + z) < -r- < 3 C (1 + z) 2 + 2 D (1 + 2) , (12) 

uncertainties in the cosmological parameters allow for both types of DE. 

One should realize that the condition (JZJ) is local: Dark Energy is of the Phantom 
type for any z where this condition is satisfied. Realization of (jJJ) even for one redshift 
value z would constitute in itself a crucial result as it would discard all models that 
can not be of that type like quintessence models (minimally coupled scalar fields), 
and this would force us to consider alternative Phantom DE models. On the other 
hand, the usefulness in this respect of the inequalities ((Zj) and (|10JllJ12j) depends 
on the accuracy with which we can recover the quantity h(z) from the data as a 
continuous function of redshift (see [T7] for a recent discussion). Indeed, the use 
of these inequalities imply diffentiating twice the luminosity distance di{z). A first 
differentiation will yield h(z) 

where a prime denotes differentiation with respect to z and we have introduced the 
dimensionless Hubblefree luminosity distance D L (z) = H d^z). The equation of 
state parameter w(z) requires a second differentiation 

w[z) v - n mfi (1 + z f - n kfi {i + z y [ > 

We recognize from eq. (|14|) the expression for the DE pressure Pde{ z ) that can be 
recovered from the data 

Pde(z) = (SnG)- 1 ((1 + Z )^--3H 2 + Q kfi H 2 (l + z) 2 ^ . (15) 
We see that wq = w(z = 0) is found from 

1+2 dh 2 I 1 1 1 

W = -^-^ 2 j (16) 

where Qde,o can be obtained from f2 mj0 ^nd f^o, namely 

Qde,o — 1 — ^m,0 — ^fc,o • (17) 

Equations (jl3H16j) generalize the corresponding flat space expressions, they apply to 
arbitrary equations of state and allow for a determination of w(z) from the di{z) 
data. 



4 



Figure 1: Several iso-di curves, curves corresponding to the same value of the lu- 
minosity distance di, are shown in the (w,Qk,o) parameter plane for arbitrary but 
constant Dark Energy equations of state. The uncertainty in the universe curvature 
is essentially in the range —0.15 < Qk,o — 0.05. It is seen that uncertainties on the 
value of Qk,o can induce significant uncertainties on the quantity w. This degener- 
acy is particularly interesting in the neighbourhood of the point (w = — l,{lk,o = 0) 
because w < — 1 (Phantom Dark Energy), w = — 1 (cosmological constant A), and 
w > — 1 all have dramatically different implications on the nature of Dark Energy. 

3 Behaviour in the parameter space (w, flk,o) 

We want to investigate the amount of cosmic confusion in the region around w = —1 
that can arise due to a non vanishing curvature of the universe. As said in the 
Introduction, we will work in the two dimensional parameter space (w,Qk,o) f° r fixed 
fl m fl and arbitrary constant w. The observational bounds on the geometry of the 
universe coming from the WMAP (Wilkinson Microwave Anisotropy Probe) data 
|10j . possibly combined with other data, are the most stringent ones. They favour a 
marginally closed, nevertheless close to flatness, universe which is reassuring for the 
inflationary paradigm. The bounds extracted from the data depend on the various 
priors about the underlying model universe and most crucially on the nature of DE. 
We certainly don't want to restrict ourselves to w = — 1 and uncertainties in the DE 
equation of state will inevitably relax the observational bounds. We will take the 
following interval of interest jTTj 



-0.15 < Vl kfi < 0.05 . 




On the other hand, the data allow w = — 1 as a good fit to the SNIa data. So it is 
clear that the interesting region is centered around (w = —l,Qk,o — 0) and we can 
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Figure 2: Curves are displayed in the (w, z cr ) plane that give for arbitrary but constant 
w the corresponding critical redshift value z cr for which 9 g£^ (z = z cr , Cl k>0 = 0) = 0. 
Each curve corresponds to a fixed value of the cosmological parameter Q m ,o and we 
have from top to bottom f2 m0 — 0.5,0.4,0.3,0.2,0.1. For the value w = —1, the 
coresponding redshift interval is 2.8 < z cr < 3.7 with z cr = 3.1 for the fiducial value 
^m,o — 0-3. At z — z cr , the degeneracy is maximal in the flk,o direction, but minimal 
in w. 



expect future data will continue to single out this region . 

We consider first luminosity-distances d^ in our two dimensional parameter space, 
measured at fixed redshift z. It is clear from the iso-c^ curves, curves having the 
same d^ value, displayed on Figure 1 for fixed z = 1.5 how cosmic confusion can 
arise: observations interpreted as a pure cosmological constant A assuming a flat 
universe can equally well be interpreted as either DE with w > — 1, resp. Phantom 
DE (w < —1), if one takes a closed, resp. open, universe. 

The crucial point is how this evolves with increasing redshift. In the region of 
parameter space of interest, we find that this behaviour changes at some critical 
redshift z cr ~ 3. Indeed at z — z cr , there is essentially no degeneracy in w while we 
have a maximal degeneracy with respect to curvature 

dd L {z = z cr ) 

dn kfl k >°=° - • (19) 

The value of z cr depends of course both on the values of fl m ,o and w. We have 
checked that a corresponding behaviour takes place when we consider different values 
0.1 < fi m ,o < 0.5 and different w, as is summarized in Figure 2 where we can read 
the critical redshift that corresponds to various values fi mi o and w. So, for given w 
of interest, our analysis can be repeated for any Q m0 as its value will be determined 
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with increasing accuracy. 

The degeneracy in the (w,Qk,o) parameter plane is opposite for higher redshifts 
z > z cr : the iso-di curve through the point (w = —1,0^0 = 0) fits also a closed 
universe with Phantom DE as well as an open universe with w > — 1. One could use 
this property in order to constrain the degeneracy in w as can be seen from Figures 
3, 4. It could also be used in order to constrain substantially the geometry of our 
universe. However one should remember that our analysis is done for non-varying 
equations of state and that the precision with which we can recover the luminosity 
distance at z > z cr is a crucial issue that cannot yet be settled. On the other hand, 
results found here will essentially hold if Dark Energy has a weakly varying equation 
of state, and they would become particularly interesting if w ~ — 1. 
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Figure 3: Two iso-d L curves (solid) are shown for f2 mj0 = 0.3 in the (u>,f2 fej0 ) plane 
that pass through the point (w = — l,fi fc0 = 0), the first (upper solid curve on the 
left) correspond to z = 1 while the second solid curve corresponds to z — 4. The 
corresponding 1% errors (dashed lines) are also shown for each curve. It is seen that 
if measurements at both redshifts fit the couple (w = —l,Qk,o = 0), the degeneracy 
in the value of w that comes from the possibility to have non-flat universes could be 
substantially restricted as well as the spatial curvature. 



4 High z measurements using gravitational waves 

High redshifts measurements of the luminosity distance di could be made using the 
emission of gravitational waves (GW) by inspiralling black holes binaries at high 
redshifts z [TH|. The flux F, the amount of energy per unit of time and per unit of 
area in gravitational waves produced by such binaries that we measure, is given by 
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(we put back the velocity of light c) 



c 



,8 



(K + K), (20) 



16ttG 



where h x are the two polarization states. Formally we can express dL, the energy 
emitted by the system per unit of time (in the binary's restframe) in the solid angle 
dQ, in the following way 



dL = T ^(hl + hl)dldn, (21) 

so that, with a the present value of the scale factor and r\ the coordinate distance 
to the binary, the measured flux F becomes 



dL dL 1 dL 



(2? 



dA(l + z) 2 a 2 Q r\ dVt (1 + z) 2 d\ dn ' 
and (f2T)J) is recovered. The polarisation states h + , h x are given by 

5 5 

h + = 2 -^~ ^ (7r/(t))I [l + (L.n) 2 } cos2$(t) (23) 

hx = ^ (7r/(t))§ (2 ^-^ sin2$(t) • (24) 

The unit vector n defines the binary's position on the sky while L points along the 
binary's orbital angular momentum and hence defines its orientation. The important 
quantity Ai is the chirp mass given by 

^=(1 + ,)-^=^, (25) 
(mi + m 2 ) s 

where m 12 are the two masses. As can be seen from ([25)1 . M. contains the factor 1+z. 
The data allow us to determine M. accurately but without independent knowledge of 
the masses m 12 , we cannot deduce the redshift z from the data. Finally the orbital 
phase $(£) is related to the wave frequency fit), ^ = n f(t) while f(t) depends on 
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the chirp mass A4, f ~ / a M.^. Extracting /, or A4, from the data allows us to 
find dt-, for example we can use the ratio j 

where g x ,+ encodes the dependence on the binary's orbital phase, position and incli- 
nation and can be determined from the data. To summarize, we can extract from the 
GW data the quantity d^, however not the redshift z that corresponds to it. 

The spatial interferometer LISA for the detection of gravitational waves around 
frequencies v ~ (1CT 3 — 10~ 4 ) Hz that will operate in the future could possibly 
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make such high z measurements. For example the gravitational waves emitted by 
inspiralling binary black holes with large masses ~ 10 5 M at redshifts z > 
z cr = 3 — 4, and even at higher redshifts, could be detected by LISA and this detection 
would yield a measurement of the distance cIl of the source but not of its redshift. If 
in addition, some electromagnetic counterpart to the GW emission can be identified, 
then the situation changes dramatically. By determining the position of the source, 
the relative error on the luminosity distance could even drop below the 1% level, 
^ ~ 1% PHI- Note that these black hole binaries could constitute standard candles 
( "sirens" ) complementary to the SNIa. 

In this way, new di points at high redshifts could be obtained with high accu- 
racy. Clearly this could considerably restrict the possible degeneracy in the (w, flk,o) 
plane and it is this aspect we want to stress here. We will assume for reference that 
measurements of di(z) at z ~ 4 are made at the 1% precision level. 



-0.8 




Figure 4: As in Figure 3, two iso-c^ curves (solid) are shown, but this time they 
fit at both redshifts z = 1 (upper solid curve on the left) and z = 4 the values 
(w = —1.1,^0 = 0). With 1% errors at both redshifts, it is seen that the equation 
of state parameter w is constrained to satisfy w < —1. 

Let us assume for example that high redshifts, as well as low redshifts luminosity 
distance measurements fit a universe with w = — 1 and f^o = 0. Then, as can be 
seen from Figure 3, the degeneracy in w would become minimal, essentially that of 
measurements made at z — z cr . It would be natural to interpret this as support for a 
flat A dominated universe. A similar reasoning would apply if observations, assuming 
a flat universe, fit a constant w different from —1: if di measurements at high and 
low redshifts are obtained that fit a flat universe with Phantom DE, the uncertainty 
on the equation of state could be small enough so as to provide strong evidence for 
Phantom DE, again assuming a constant w, as we show with an example in Figure 
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4. We see also from Figure 4 that we must have w < —1.06 in order to characterize 
phantom Dark Energy unambiguously. 

It would be interesting to investigate how a varying equation of state would mod- 
ify these conclusions and we leave this for further investigation. Clearly, such an 
investigation will become more effective once we have enough information about the 
time- varying equation of state that could be extrapolated to redshifts z > z cr . It is 
clear however that our results will essentially hold for a weakly varying equation of 
state. In this case accurate results can only be obtained for given equation of state. In 
conclusion, we have shown how the complementarity of low redshift SNIa data with 
high redshift measurements by LISA could provide us with useful information on the 
nature of Dark Energy. We feel this complementarity is an interesting direction to 
explore in particular if one relaxes the assumption of spatial flatness. 
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